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ABSTRACT 

We present a model for pulsar radio eclipses in some binary systems, and test this model for 
PSRs B1957 + 20 and J2051 — 0827. We suggest that in these binaries the companion stars are 
degenerate dwarfs with strong surface magnetic fields. The magnetospheres of these stars are 
permanently infused by the relativistic particles of the pulsar wind. We argue that the radio waves 
emitted by the pulsar split into the eigenmodes of the electron-positron plasma as they enter the 
companion's magnetosphere and are then strongly damped due to cyclotron resonance with the 
ambient plasma particles. Our model explains in a natural way the anomalous duration and 
behavior of radio eclipses observed in such systems. In particular, it provides stable, continuous, 
and frequency-dependent eclipses, in agreement with the observations. We predict a significant 
variation of linear polarization both at eclipse ingress and egress. In this paper we also suggest 
several possible mechanisms of generation of the optical and X-ray emission observed from these 
binary systems. 

Subject headings: binaries: eclipsing — pulsars: individual (PSR B1957 -I- 20, J2051 — 0827) 



1. Introduction 

Eclipsing millisecond pulsars were expected to 
be a missing link in the evolutionary connec- 
tion between millisecond pulsars and low-mass X- 
ray binaries (Alpar et al. 1982). To the best of 
our knowledge, the following eclipsing millisec- 
ond pulsars have been observed for today: PSR 
B1957-F20 (Fruchter et al. 1988), PSR B1744-24A 
(Lyne et al. 1990), PSR B1259-63 (Johnston ct al. 
1992), PSR J2051 - 0827 (Stappers et al. 1996a), 
47 Tuc I, J (Manchester et al. 1991), O, and R 
(Camilo et al. 1999). Up to now, PSRs B1957 + 20 
and J2051 — 0827 are the best studied of them. 
These two pulsars exhibit stable, anomalously long 
and frequency-dependent radio eclipses, that is, 
eclipses are long at low observed frequencies, and 
become shorter at higher frequencies. 

In this paper we propose a physical mechanism 
for the eclipsing binary pulsars and test it for 
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PSRs B1957 + 20 and J2051 - 0827. Our model 
assumes that the companion stars are magnetic 
white dwarfs. Relativistic particles from the pul- 
sar wind are trapped by the companion's mag- 
netic field. This leads to formation of the ex- 
tended magnetosphere of the companion star. We 
demonstrate that the pulsar radio emission un- 
dergoes strong cyclotron damping while passing 
through the companion's magnetospheric plasma. 
The model accounts for the most peculiarities of 
the observed eclipse pattern in both aforemen- 
tioned eclipsing systems. To begin with, let us 
shortly review the main observational features of 
these interesting objects. 

1.1. PSRB1957-h20 

One of the two fastest pulsars known (with 
the period P « 1.6 ms and the spin-down rate 
P « 1.7 X 10^2°), the "black widow" pulsar 
B1957 + 20 shows regular and entirely periodic 
eclipses at meter wavelengths, which occupy about 
10% of its Torb « 9.2 hr orbital period. The 
eclipses are quite stable in length at any given ob- 
serving frequency, although their length depends 
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strongly on the frequency: at 318 MHz it aver- 
ages to about 55 minutes, but decreases to about 
33 minutes at 1.4 GHz. This frequency depen- 
dence appears to fit a Ate oc ;y-0-4±o.i pQ-^^Qj- 
law, where Ate denotes eclipse duration at a fre- 
quency u. The orbit of the binary system is rather 
circular, with the separation of the companions 
a « 1.7 X 10" cm « 2.5 i?© (Brookshaw & Ta- 
vani 1995), where Rq = 6.96 x 10^° cm is a solar 
radius. 

At lowest observing frequencies eclipse ingress 
is quite rapid, whereas eclipse egress is slower and 
somewhat turbulent. At 1.4 GHz, however, this 
asymmetry is not observed. The excess electron 
density on either side of the eclipse is found to vary 
by a factor of two from eclipse to eclipse. The de- 
lay between the left and right circularly polarized 
pulse arrival times were also detected (Thorsett 
et al. 1989; Pruchter et al. 1990; Ryba & Taylor 
1991). 

HST optical images of the PSR B1957-F20 made 
at different orbital phases display a dramatically 
variable system which is brightest when the side of 
the companion heated by the pulsar wind faces the 
Earth, and darkens by several magnitudes when 
the opposite cool side comes into view (Fruchter 
1995). Companion's radius, i?opt ~ 0.12 i?©, was 
derived from the optical images using the esti- 
mated distance to this object d = 0.8 kpc. There- 
fore i?opt is less than one-half the Roche lobe ra- 
dius Rl « 0.29 i?0, and the enigmatic companion 
of PSR B1957-h20 does not even fill its Roche lobe 
(Djorgovski & Evans 1988). 

Soft X-rays from this system were detected 
using the ROSAT observatory (Fruchter & Goss 
1992; Kulkarni et al. 1992). An X-ray luminos- 
ity of about 10^^ erg (which is ~ 10~^ of the 
pulsar spin-down energy) was found, although no 
variability has been detected. Below we discuss 
possible sources of this high-energy and optical ra- 
diation. 

1.2. PSR J2051 - 0827 

The discovery of this eclipsing millisecond bi- 
nary system was reported by Stappers ct al. 
(1996a) who carried out timing observations at 
several frequencies between 408 MHz and 2.0 GHz. 
The orbital period Torb ~ 2.38 hr of this system 
makes it the third shortest pulsar orbital period 



known, behind 47 Tuc R (Camilo et al. 1999) 
and PSR B1744-24A (Lyne et al. 1990). Thus, 
the system is extremely compact, with a separa- 
tion of the binary components a « 1.0 i?Q. High- 
precision timing observations of PSR J2051 — 0827 
(the dynamical parameters of this pulsar are: 
P = 4.5 ms and P = 1.3 x 10"^°) indicate 
that the orbital period is decreasing at a rate of 
To,b = (-11 ± 1) X 10"^^ (Stappers et al. 1998). 
Such an orbital period derivative implies a decay 
time for the orbit of only 25 Myr, which is much 
shorter than the expected timescale for the abla- 
tion of the companion. This, in combination with 
a very slow transverse velocity of the pulsar it- 
self, questions the formation of the system in the 
standard manner (Stappers et al. 1998). 

The duration of the pulsar eclipse at 436 MHz 
is ^ 10% of the orbital period, which implies a ra- 
dius of the eclipse region surrounding the compan- 
ion i?e « 2 X 10^° cm « 0.3 i?0. The variation in 
eclipse duration between 436 and 660 MHz seems 
to be frequency dependent with Aig oc j/""--^^. On 
the other hand, at 1.4 GHz the pulsar emission 
was detected throughout the low- frequency eclipse 
region in approximately half of the observing ses- 
sions (Stappers et al. 1996a). These results, along 
with the observations at 2.0 GHz, suggests that no 
eclipse occurs at high radio frequencies. 

Stappers et al. (1996b) also detected an opti- 
cal counterpart of the companion of PSR J2051 — 
0827. The amplitude of its light curve is at least 
1.2 mag. Thus, PSR J2051 - 0827 is the sec- 
ond pulsar binary for which direct heating of the 
companion was observed. The radius of the com- 
panion's optical counterpart (derived in the black- 
body approximation) falls in the range i?opt ~ 
0.067 - 0.18 i?Q (Stappers et al. 1996b). Thus, 
the companion star is likely to fill its Roche lobe 
(i?i « 0.1'SRq). The best-fit models, based on 
the new photometry of the secondary star (Stap- 
pers et al. 1999), require that greater than 30% of 
the incident energy is absorbed by the companion 
and reradiated as optical emission. The unillu- 
minated side of the companion (i.e., observed at 
an orbital phase corresponding to the pulsar ra- 
dio ecfipse) was found to be cool (Stappers et al. 
1999), with a best-fit temperature of 3000 K, sim- 
ilar to that obtained for the companion to PSR 
1957-F20 (Pruchter et al. 1995). 
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2. Nature of the companion stars 

Let us begin the discussion with a few com- 
ments on the nature of the companion stars. The 
minimum masses of the companions in both cases 
were calculated from the mass function /{irip, rric) 
using observed orbital parameters and assum- 
ing the pulsar mass nip « 1.4 and the or- 
bital inclination i > 60°, the latter necessary 
for the existence of eclipses. For the compan- 
ions of PSRs J2051 - 0827 and B1957 + 20 this 
yields nic ~ 0.027 Af© (Stappers et al. 1996a) 
and rUc « 0.022 Mg (Fruchter et al. 1988), re- 
spectively (i.e., very low companion masses). The 
actual masses of companions are within a factor 
of 1.2 of the above minimum values (because of 
i > 60°). Fruchter (1995) argues that the com- 
panion of PSR 1957-1-20 is obviously well below 
the hydrogen burning limit, and one would expect 
this star to be a degenerate dwarf. We suggest 
that the same is valid for the companion of PSR 
J2051 - 0827. The radius of a low-density white 
dwarf (assuming nonrelativistic electron Fermi gas 
in its interior) is given by the mass-radius relation 
(Shapiro & Teukolsky 1983) 

^.4.05.10-(^)""'%-'3, (,) 

where /Xe = A/Z is a mean molecular weight per 
electron^. Substituting the masses of each com- 
panion star in equation (1) we can calculate their 
radii for the cases of pure hydrogen {p,,, = 1) 
and pure helium (/Xe = 2) white dwarfs, which 
yield Rf^ = 0.145 i?©, Rf^ = 0.046 i?©, R^2 = 
0.135 i?0, and Rf^ = 0.043 i?©, respectively. For 
a mixture of 25% He and 75% H by mass {X = 
0.25), we obtain that i?ci ~ Rc2 ~ 0.1 Rq. We see 
that these values of companion radii arc compa- 
rable to those estimated from the optical observa- 
tions (see Section 1). 

Because magnetic flux is conserved in gravi- 
tational collapse, the stellar magnetic fields can 
be amplified to more than a million Gauss when 
a star contracts to the white dwarf stage. Ob- 
servations show that the magnetic fields at white 
dwarf surfaces vary from a few tens to millions of 

^For brevity, hereafter subscripts "1" and "2" correspond to 
the physical quantities regarding the companions of PSRs 
B1957 -I- 20 and PSR J2051 - 0827, respectively. 



Gauss (Lang 1992). We assume that both white 
dwarfs discussed above possess significant mag- 
netic fields at their surfaces, and estimate these 
magnetic fields below. 

3. Eclipse mechanism 
3.1. Review 

Thompson et al. (1994) and Thompson (1995) 
analyzed a variety of physicial mechanisms which 
can cause pulsar eclipses. They considered eclipses 
by a wind from the stellar companion, by a stellar 
magnetosphere, or by material entrained in the 
pulsar wind. The refractive eclipse proposed by 
Phinney (1988) was ruled out by Thompson et al. 
(1994) due to inability of this mechanism to ex- 
plain a sensitive dependence of the eclipse dura- 
tion on frequency, as well as measured time de- 
lays. The other mechanisms discussed by Thomp- 
son et al. (1994) and Thompson (1995) are as fol- 
lows: a free-free absorption, scattering by plasma 
turbulence, induced Gompton scattering, stimu- 
lated Raman scattering parametric instability and 
synchro-cyclotron absorption. Thompson et al. 
(1994) and Thompson (1995) are inclined to the 
opinion that the eclipse is due to cyclotron absorp- 
tion by plasma of temperature Te ~ (1 — 4) x 10^ 
K and density rig — 4 x 10^ cm~^ in a field of 
strength Be ~ 15 — 20 G. Another possibility is a 
synchrotron absorption by relativistic plasma elec- 
trons (Eichler 1991). 

Detailed discussion of eclipse mechanisms was 
also presented by Eichler & Gedalin (1995). They 
focus on three wave processes, in particular Ra- 
man scattering (the decay photon — > photon + 
plasmon stimulated by an ambient plasmon field) 
and induced plasmon emission (i.e. the same pro- 
cess stimulated by photons already in the final 
state). These authors claim that both mechanisms 
are able to produce pulsar eclipses via large an- 
gle scattering on the initially beamed pulsar radi- 
ation, and neither of them requires high plasma 
densities. Brookshaw & Tavani (1995) considered 
different scenarios of the mass outflow in bina- 
ries, such as mass outflow intrinsic to the com- 
panion star and pulsar irradiation-driven outflows 
or Roche lobe overflow, and discarded the latter. 

A mechanism based on three-wave interactions 
involving low-frequency acoustic turbulence was 
proposed to be responsible for pulsar eclipses by 
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Luo & Melrose (1995); Luo & Chian (1997). The 
authors claim that the eclipsing material consists 
of plasmas from the companion wind, heated by 
the pulsar wind. This may result in a hot corona- 
like plasma which is possibly non-isothermal. 
They considered both the case of electron acous- 
tic waves when the ion acoustic temperature Ti is 
higher than the electron temperature Tg and that 
of ion acoustic waves when T^, > Ti. It was found 
by Luo & Chian (1997) that induced scattering 
off electron acoustic waves can be important and 
can cause pulsar eclipse. At the same time, it was 
demonstrated that induced Brillouin scattering 
involving low-frequency ion acoustic waves can- 
not be the main reason of eclipse, because of the 
relatively small growth rates of these waves. 

3.2. Geometry of the model 

As we have seen in Section 1, the ratio of the 
eclipse duration (at low frequenices) to the orbital 
period is about the same in both PSR B1957-|~ 20 
and PSR J2051 - 0827, i.e., Aie/Torb ~ 0.1. As 
the eclipse nears, the pulse decreases in ampli- 
tude but does not change in shape, showing that 
the eclipse mechanism removes photons from the 
line of sight but does not significantly scatter the 
pulse. Thus, only absorption or large-angle scat- 
tering may be responsible for it (Stappers et al. 
1996a). On the other hand, the eclipse length 
indicates that the eclipsing medium extends to 
the distances much larger than the dimensions 
of a companion of any conceivable composition 
(Fruchter 1995). Indeed, for PSR B1957-H 20 the 
radius of eclipsing region i?e(318MIIz) w 0.76 i?© 
(so that the "eclipsing spot" hes at least 1.5 i?© 
across), i?e(l-4GHz) « 0.46 i?©, whereas for PSR 
J2051-0827i?e(436MHz) «0.31i?©. Weseethat 
in both cases (at all radio frequencies) the eclipse 
region lies even well beyond the corresponding 
companions' Roche lobes {Rli ~ 0.29 i?© and 
Rl2 ~ 0.13 i?©, respectively). 

Thompson (1995) points out that the large 
size of the eclipses (compared to the companion's 
Roche lobe) could be explained either by diffusion 
of plasma into the pulsar wind, or by a large com- 
panion magnetosphere. We adopt this statement 
and claim that companions of both pulsars are 
degenerate dwarfs with significant magnetic fields 
at their surfaces. At the same time, we suggest 
that the magnetospheres of these magnetic white 



dwarfs are infused with relativistic particles sup- 
plied permanently by the pulsar relativistic wind. 
Figure 1 represents schematically geometry of such 
a binary system. 




Fig. 1. — Schematic representation of a binary sys- 
tem (viewed at some angle) at an orbital phase 
corresponding to the radio eclipse, a is the bi- 
nary separation. An extended magnetosphere of 
the white dwarf is continuously supplied by rela- 
tivistic plasma particles of the pulsar wind. 

Let us estimate the density of a wind plasma at 
the distances from the pulsar surface correspond- 
ing to the binary separation a. The number den- 
sity of pair plasma inside the pulsar light cylinder 
(with the radius i?Lc = Pc/2'k) is 
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where k is a Sturrock multiplication factor (Stur- 
rock 1971) and a is an inclination angle between 
the pulsar magnetic and rotation axes. From equa- 
tion (2) we can calculate the number density of 
particles at distance of one light cylinder radius, 
rip^ . Assuming further that the density falls ac- 
cording to the inverse square law beyond the pul- 
sar light cyhnder n^{R) — n^^{R^^^/ R)"^ we get 
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(3) 

Substituting in equation (3) the dynamic param- 
eters of each of the pulsars discussed in this pa- 
per, taking a ^ 90° for PSR B1957 + 20 (note 
that this pulsar has an interpulse indicating that 
it is an almost orthogonal rotator), and assum- 
ing that the inclination angle is rather large for 
PSR J2051 - 0827 too, we find that at the dis- 
tances corresponding to binary separations in each 
of the systems (ai « 2.5 i?© and 02 ~ 1.0 i?©, re- 
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spcctivcly), the plasma densities in pulsar winds 
are n^{ai) « 3.1 x 10~^k cm~^ and n^{a2) « 
3.6 X 10~^K cm~^, respectively. We see that the 
densities are about the same order of magnitude 
in these two cases. Note that the actual value 
of plasma density in the companions' magneto- 
spheres may be much higher than inferred from 
equation (3). The point is that the slowest par- 
ticles of pulsar wind trapped by the companion's 
magnetic field should keep bouncing between its 
magnetic poles without relative density loss on 
an extended timescale. Such an accumulation 
of plasma particles may apparently increase their 
number density. 

3.3. The model 

We share the opinion expressed by Thompson 
et al. (1994) that the most promising model of 
eclipsing pulsar binaries is cyclotron absorption. 
However, we believe that the proper considera- 
tion of waves damping in the relativistic magne- 
tized plasma should be performed basing on ki- 
netic treatment of the cyclotron damping. Bas- 
ing on this statement, as well as on the physical 
picture of an eclipsing binary system presented in 
Section 3.2 (sec also Fig. 1), we presented a model 
for the radio eclipse of PSR B1957 + 20 (Khechi- 
nashvili & Melikidze 1997). Here we develop and 
further extend this model by applying it to PSR 
J2051 - 0827. 

Let us check the possibility of waves damping 
in the plasma of companion's magnetosphere. Ra- 
dio waves emitted by the pulsar (with the vacuum 
spectrum u> = kc) should split into the eigenmodes 
of the relativistic pair plasma, as they enter the 
companion's magnetosphere. It was shown (e.g., 
Lominadze et al. 1986; Arons & Barnard 1986) 
that in the magnetized relativistic pair plasma 
there exist three wavemodes with the frequencies 
much less than ojj^/'j^, where ui^ = eB/mc is the 
gyrofrcquency and 7 is a mean Lorentz factor of 
the plasma particles. In the general case of oblique 
propagation with respect to the local magnetic 
field, these eigenmodes are: i) the purely electro- 
magnetic extraordinary (X) mode; ii) the sublumi- 
nous Alfven (A) mode; and iii) the superluminous 
ordinary (O) mode. The last two modes are of 
mixed electrostatic-electromagnetic nature. Elec- 
tric field vectors E*-* and of O and A-modes 
lie in the (k B) plane, while the electric field of X- 



mode is directed perpendicularly to this plane. 

All the three wavemodes may be strongly 
damped in the companion's magnetosphere at the 
cyclotron resonance 



bj - ki\v = Q 

/ res 



(4) 



on the particles of ambient plasma with the mean 
Lorentz factor 7^^^ = 7^ . Note that in equation (4) 
k\\ denotes the component of the wavevector along 
the local magnetic field, and we neglected the cur- 
vature drift of the bulk plasma particles across the 
plane of the magnetic field line curvature. 

Let us calculate the characteristic frequency 
of damping of the plasma eigenmodes at the cy- 
clotron resonance. We assume that all the modes 
are vacuum-like in the domain of their spectrum 
where the cyclotron damping occurs, i.e. their 
spectrum can be represented as 



u) = kc{l — 6), where (5^1. 



(5) 



Using the dispersion relation in this form we can 
rewrite the resonance condition (4) as follows 



kc{l-6) -kccos0 il- ^] ='^, (6) 
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where is an angle between the wave-vector and 
the magnetic field in the resonance region, and we 
used V Ki 1 — 1/27^. Prom equation (6) we obtain 
two approximations of the frequency of damped 
waves: 

u;,^2-f^to^, if e^ — , (7) 
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and 



7^(1 -cos (9) 7^ 

The decrements of cyclotron damping for the 
plasma eigenmodes are obtained by solving the 
corresponding dispersion equations (see, e.g., 
Lominadze et al. 1986; Khechinashvili 1999) for 
the imaginary part of the complex frequency 
{u! ^ oj + ir). Therefore the decrements of X, 
O and A-modes write 



(9) 
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respectively, where 
g (6) = cos^ e + 
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and 7y is a thermal spread of the particles energy. 
In equations (9) and (10) we used the mean value 
theorem and the normalization of the distribution 
function, namely that / (7^) 7r = / / (7) ^7 = 1. 
Note that for small angles (or small Lorentz factors 
of the ambient plasma particles) g (6*) « 1 and the 
decrement of O and A- modes (eq. [10]) transforms 
into that of X- waves (eq. [9]), 
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Fig. 2. — Angular dependence of the decrement of 
O and A- modes in the limit j 9 ^ 1. The function 
g{9) defined by equation (11) is plotted versus 9 
in degrees, for three different values of the mean 
Lorentz factor - 7^^ = 2, 7^ = 10 and 7^ — 100 - 
represented by the dotted, dashed and solid lines, 
respectively. 

In the opposite limit, i.e. when the angles (or 
the plasma mean Lorentz factors) are relatively 
large, g{6) is plotted in Figure 2, for different val- 
ues of 7p . It is seen that if 1 /^^ <^ 9 < 1 the value 
of the decrement of O and A-modes is a few times 
larger than that of the X-mode (provided that the 
other parameters remain the same in both cases). 
For example, for 9 = 30° and 7^ — 100 (a dashed 
line in Fig. 2) we have r(o,A)/rx ~ 4. The two 



decrements become equal each other {g{9) « 1) for 
the angle 9 ~ 70°. On top of that, it appears that 
O and A-modes propagating nearly transversely 
to the local magnetic field undergo no cyclotron 
damping, since g (0) ^ if — > 90°. Indeed, the 
A-mode vanishes at exactly transverse propaga- 
tion to the local magnetic field (its group velocity 
tends to zero), whereas the 0-mode transforms 
into a mode whose electric field is parallel to the 
external magnetic field and k ± E. Such a wave 
is not subject to cyclotron damping. 

It is obvious that the angles between wavevec- 
tors and the local magnetic field cannot be small, 
hence, 9 ^ l/7p- Therefore, we can calculate 
the frequency of damped waves from equation (8), 
which is now convenient to represent in the follow- 
ing form 
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Here r is a distance from the companion's center, 
Rc is the radius of the companion star, and Be is 
its surface magnetic field measured in Gauss. Note 
that we assume the companion's magnetic field to 
be dipolar, i.e. B (r) = Bc{Rc/r)^- 

V >v' 



/ \ v>v" 


r 




V >V" 

^ 



V >v 



Fig. 3. — "Eclipsing spots" at different frequen- 
cies. The pulsar emission reaching the damping 
region represents a full range of radio frequencies. 
It is seen that all the frequencies lower than v' are 
damped at distances from the stellar surface larger 
than r' . The waves with higher frequency v" > v' 
propagate freely through this outer region with a 
low value of the magnetic field strength B' . How- 
ever, they are damped at the distance r" where 
magnetic field value is higher, B" > B' . All the 
radio waves with frequencies v > u" escape this 
narrow inner region. 
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It is seen from equation (13) that the frequency 
of damped waves is inversely proportional to the 
cube of the distance from the companion's sur- 
face (obviously, this dependence results from i^^ oc 
U!g^ oc Be, the latter being the magnetic field value 
in the eclipse region). In other words, distinct fre- 
quencies are damped at distinct heights above the 
stellar surface (corresponding to appropriate val- 
ues of the magnetic field), while higher frequen- 
cies are damped closer to the companion star. The 
waves with frequencies somewhat higher and lower 
than the one given by equation (13) propagate al- 
most freely in the plasma at this distance from 
the star. Although, they can also reach a re- 
gion at a different altitude or with different direc- 
tion of local magnetic field (hence, angle 9) and 
be damped there. Figure 3 represents schemati- 
cally a top view of the eclipse region. It follows 
from this picture that the size of "eclipsing spot" 
is frequency-dependent, being larger at lower fre- 
quencies. Indeed, radio waves with a low fre- 
quency i>' (from the entire range of the pulsar radio 
spectrum) are damped at some large distance r' 
from the white dwarf surface. At the same time, 
waves with higher frequency ly" > z/' propagate 
almost freely through this outer region, although 
they are damped as they reach the altitude r" < r' 
from the stellar surface with a higher value of mag- 
netic field. 

Let us find the values of physical parameters 
matching the observational data on PSR 1957-1-20 
with our eclipse model. Substituting first = 
0.318 GHz and then i/^ = 1.4 GHz in equa- 
tion (13), and using corresponding values of eclipse 
region radii at these frequencies r = Rei^j) we can 
find that this range of radio frequencies is damped 
at the distances 0.46 — 0.76 Rq from the stellar sur- 
face, if 



,f ^ 50 R%. 

7p(l-cos6») 



(14) 



First assume that the companion star is a hydro- 
gen white dwarf with Rc = 0.145 J?©. Speculating 
further that the mean Lorentz factor of the par- 
ticles of companion's magnetospheric plasma are 
moderate 7^ ~ 10 — 100 and the average angle is 
rather large, 6 = 60° — 90°, we obtain that the 
magnetic field strength at white dwarf's surface 
Be 2 X 10^ - 3 X 10^ G. Alternatively, if the 
companion is built by pure helium (so that Rc = 
0.046 Rq), equation (14) yields ~ 5 x 10<5 - 10^ 



G for the same set of parameters. 

Applying the above procedure to PSR J2051 — 
0827 system we obtain that the band of radio fre- 
quencies up to about 2 GHz is damped in the com- 
panion's magnetosphere, if the following condition 
should is satisfied. 
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Here it was taken into account that the waves with 
= 436 MHz are damped at the distance Rc = 
0.31 i?0. Assuming, as in the previous case, that 
7p ~ 10 - 100 and 6* = 60° - 90°, we can estimate 
that Be ~ 2x10^^-4x10^ G for the hydrogen white 
dwarf (i?c = 0.135 iio), and B,, - 5 x 10^ - 10^ G 
for the helium white dwarf {Rc = 0.043 i?©). 

For all the intermediate cases of white dwarfs 
made of various mixtures of He and H by mass 
(that is, various values of X, see Section 2) the sur- 
face magnetic field falls in the range between the 
above values of Be in both cases. Detailed spectral 
observations of both companion stars would deter- 
mine the spectral class (hence, the chemical com- 
position) of each white dwarf, as well as constrain 
the values of surface magnetic fields at their sur- 
faces. In the frame of our model (eqs. [14,15]) this 
would provide an indirect measurement of average 
angles 9 and mean Lorentz factors of companion's 
magnetospheric plasma particles. The latter, in 
turn, may specify the energy of the slowest (and 
the most numerous) particles of the pulsar wind. 
All in all, we see that both companion stars pos- 
sess quite strong surface magnetic fields, hence we 
deal with so called magnetic white dwarfs, proba- 
bly with mega-Gauss magnetic fields. Let us no- 
tice that such objects are quite common in the 
Galaxy (Lang 1992). 

The necessary condition of the cyclotron res- 
onance development is that the characteristic 
timescale of the instability, = l/lFdj, should 
be much less than the time waves take to escape 
the resonance region ~ 2Rc/c. The decre- 
ment Fd is given by equation (9) for X-mode, and 
equation (10) for O and A-modes. Assuming a 
relatively broad distribution function for compan- 
ion's magnetospheric plasma (so that 7^ ~ 7^) 
the former can be rewritten as follows 



= -O.i 



[GHz 



(16) 
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while the latter differs from equation (16) by an 
angle-dependent multiplier g{0) (eq. [11], see also 
Fig. 2). Thus, for pulsar to be invisible at some 
frequency h'^(Re), the following two conditions 
should be satisfied at the corresponding distance 
i?e from the companion's surface 



Te _ 2|rrf|i?e 



»1, 



(17) 




2 ir^l Re 



(18) 



Combining equations (16) and (3), expressing 
i?e from equation (13) by use of first equation (14) 
and then equation (15), we can rewrite the condi- 
tion for X-mode damping (eq. [17]) as 



6 X 10" 



% V[GHz] 



-4/3 



> 1 



(19) 



in the magnetosphere of the companion to PSR 
B1957-I-20, and as 



3 X lO"'' — 



[GHz 



-4/3 



> 1 



(20) 



in that of companion to PSR J2051 - 0827, re- 
spectively. According to equation (18), the corre- 
sponding conditions for O and A modes in each 
case differ from equations (19) and (20) by the 
angle-dependent multiplier g{6) (eq. [11]) on the 
left-hand side of these inequalities. 



- ' 1 ' 1 ' 1 ' 1 ' - 

T,=5o : 


- ' 1 ' 1 ' 1 ' 1 ' - 

: 7,= 100 






- , 1 , 1 , 1 , 1 , - 





2 3 4 1 2 3 4 

vtGHz] 



Fig. 4.— Eclipse condition for PSR B1957 + 20 
(eq. [19]). See explanations in the text. 



Fig. 5.— Eclipse condition for PSR J2051 
(eq. [20]). See explanations in the text. 



0827 



Figs. 4 and 5 display damping conditions ver- 
sus frequency in both cases. The solid lines cor- 
respond damping of X-mode at 7^ = 50 and 
7p = 100, and the dashed lines represent damping 
condition of coupled modes for the same values of 
mean Lorentz factor and the angle values 9 = 60° 
and 9 = 80°. It is clear that all the intermediate 
cases fall between the curves presented in these 
figures. The dashed horizontal line corresponds 
to the characteristic value r /r^ = 3 at which the 

e / d 

intensity of waves drops by about e^'^'^l'^d k, 400 
times. 

It is seen in Figure 4 that in the case of PSR 
B1957 -|- 20 the damping conditions are well sat- 
isfied for all the plasma wavemodes if 7^ = 50 
(hence, k = 2 x 10**). So damping is very strong 
and the pulsar radio emission certainly cannot 
reach an observer for this set of plasma param- 
eters. However, the calculations corresponding 
7p = 100 (hence, k = 10"^) indicate that for such 
parameters the pulsar emission (at least one of the 
modes) may appear visible at relatively high fre- 
quencies during a low- frequency eclipse. However, 
the present lack of high-frequency eclipse observa- 
tions of this pulsar do not allow us to judge strictly 
parameters of the magnetospheric plasma of its 
companion. 

Figure 5 shows that the set of parameters 7^ = 

100, K = lO'' and 6* = 60° - 80° are certainly ruled 
out in the case of PSR J2051 — 0827, because the 
predictions are not consistent with the observa- 
tional data from this binary system. However, it 
is seen that for = 50 (hence, k = 2 x 10^) 
some critical frequency (y^ » 1.4 GHz) may exist 
at which O and A-modes propagating with large 
angle 9 ~ 80° to the magnetic field escape from 
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the eclipse region without significant loss, while X- 
mode is still strongly damped. It means that the 
pulsar remains visible at this frequency, although 
only the waves with one polarization can be de- 
tected. If such a situation appears at the edge of 
the eclipse region (where the angle 6 varies in quite 
a small range), one can observe significant linear 
polarization at eclipse ingress and egress through 
quite a large frequency range (up to about 3 GHz, 
according to Fig. 5). On the other hand, non- 
stationary processes taking place in the compan- 
ion's magnetospheric plasma (such as, e.g., conal 
or cyclotron instabilities) may induce significant 
variations of plasma density in the companion's 
magnetospheric plasma. This may affect fulfill- 
ment of damping conditions, resulting in the even- 
tual cyclotron damping of all the wavemodes at 
the critical frequency w 1.4 GHz. Thus, one 
observes an alternating eclipse at this frequency. 
Indeed, as it was mentioned in Section 1.2, PSR 
J2051 - 0827 was detected at 1.4 GHz through- 
out the low-frequency eclipse region in about half 
of observing sessions (Stappers et al. 1996a). An- 
other reason of such an absence of eclipse at fre- 
quencies > 1.4 GHz could be a misalignment of 
the binary orbital plane with an observer's line of 
sight. Indeed, even at moderate inclination angle 
it is possible that observer misses completely the 
narrow high-frequency "eclipsing spot" (with the 
radius i?e(l-4 GHz) w 0.22 i?0), while at lower 
frequencies the pulsar still remains eclipsed at the 
same orbital phases (see Figs. 1 and 3). 

4. Discussion 

We acknowledge that the model presented 
above is to some extent idealized. This regards, 
first of all, the energy distribution of plasma parti- 
cles, which we assumed to be stable and constant 
with the altitude. In fact, a more realistic ap- 
proach would imply that more energetic particles 
accumulate at larger distances from the stellar 
surface, whereas mean Lorentz factors of plasma 
particles decrease closer to the star, due to sub- 
stantial synchrotron losses in the stronger mag- 
netic field there. It is clear that this should affect 
calculation of both damping frequency (eq. [13]) 
and decrement (eq. [16]). Indeed, according to our 
model, the eclipse length scales with the frequency 
as Ate oc On the other hand, observa- 

tions give Ate oc i/-0-4±o i for PSR B1957 -|- 20 



and Ate oc i/-" !^ for PSR J2051 - 0827, respec- 
tively. Such a deviation of the eclipse pattern from 
the theoretical prediction in both cases would re- 
sult from the aforementioned variation of plasma 
distribution function with the altitudes from the 
stellar surfaces. 

Another factor which can have a strong impact 
on the eclipse process described above is magnetic 
field orientation and its configuration. One can 
speculate about the orientation of the compan- 
ion's magnetic field with respect to the star it- 
self, and its magnetic axis could be inclined to 
the rotation axis with an arbitrary yet constant 
angle (see Fig. 1). On the other hand, the long 
history of tidal interaction with the neutron star 
in such close binary systems as PSRs B1957-I- 20 
and J2051 — 0827 are, would align the angular mo- 
menta of the white dwarf and the neutron star 
with the binary orbital angular momentum. The 
same would also bring to the situation when the 
companion faces the pulsar always with one side 
(like the moon to the Earth), a fact apparently 
confirmed by the optical observations of both sys- 
tems (see the discussion below). This implies that 
the companion's magnetosphere in both cases is 
strictly oriented with respect to the pulsar wind, 
as well as its radio emission cone (the latter being 
beamed around the pulsar magnetic axis /i). This, 
according to our model, would provide relative sta- 
bility of the eclipse pattern, excluding eclipse vari- 
ations associated with the eventual precession of 
the companion's magnetic axis around its rotation 
rotaion axis. 

Regarding the configuration of companion's 
magnetic field, we assiimcd the latter being purely 
dipolar, although one can hypothesize that the 
global magnetic field of the white dwarf is modified 
by cither a substantial contribution of local mag- 
netic fields close to the stellar surface or, e.g., a 
quadrupole component. Naturally, such a change 
of the magnetic field value and direction should 
affect waves damping in the component's mag- 
netosphere, especially at high radio frequencies 
(which arc damped closer to the stellar surface). 

It is worth estimating how the pulsar wind can 
influence the companion's magnetic field structure 
in the eclipse region. For that it is necessary to 
calculate the pulsar wind magnetic pressure 

K TT I P 
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at the distance of companion a (here I « 10'*''' g 
cm^ is the neutron star moment of inertia and R 
is a distance from the pulsar), and to compare it 
with the companion's magnetic pressure 

-c-g=4xlO-B,^(|i)', (22) 

in the ccUpse region. For the parameters of each 
of the pulsars (see Section 1) equation (21) yields 
Wp(ai) w 14 erg cm~^ and Wp{a2) ~ 3 erg cm~^, 
respectively. On the other hand, at the dis- 
tances from the white dwarf's surface correspond- 
ing to low-frequency eclipse region in each case 
(i?e(318 MHz) = 0.76 i?Q and i?e(436 MHz) = 
0.31 i?0, respectively), we find that for the typ- 
ical values given by equations (14) and (15) the 
magnetic pressure Wci = 1.3 x 10^ erg cm""^ and 
Wc2 = 2.8 X 10^ erg cm~^, respectively. We see 
that Wp <C Wc in the eclipse region in both cases, 
so it can be said that the eclipse region of white 
dwarf's magnetosphere does not "feel" the pulsar 
wind at all, and the only (yet very important) role 
of the latter is to supply this magnetosphere with 
relativistic charged particles. 

Moreover, one can estimate that the compan- 
ion's magnetosphere dominates the pulsar wind 
magnetic field in the vast expanse of both bi- 
nary systems. For PSR B1957 + 20 the "inner 
front" of the companion's magnetosphere, that is, 
the point between the white dwarf and the neu- 
tron star where Wp w Wc, is located at the dis- 
tance Tinner ~ 2.3 i?0 from thc companion, i.e. 
only about 0.2 Rq from the pulsar surface (note 
that this distance is still equivalent to about 1800 
light cylinder radii of this pulsar). Thc latter 
means that the pulsar wind particles soon fall 
under control of the companion's magnetic field. 
At thc same time, companion's magnctoitail is 
stretched up to router ~ 8.3 Rq away from the 
white dwarf. In the transverse direction com- 
panion's magnetosphere takes over pulsar wind 
up to r_L « 7.5 i?0. The corresponding dimen- 
sions of PSR J2051 — 0827 companion's magneto- 
sphere yield rjnncr ~ 0.97 i?0, router ~ 5.8 i?0 and 
rj_ ~ 5.4 i?0. We see that in this binary system 
pulsar wind particles start moving in the compan- 
ion's field after they flee from thc pulsar magne- 
tosphere on the distances more than about its 10 
light cylinder radii. 

Let us finish the discussion by a few remarks 



concerning possible physical mechanisms of the 
observed high-energy emission from these remark- 
able astrophysical systems. As it was mentioned in 
Section 1, highly variable optical counterparts of 
the companion stars were observed in both cases, 
and X-rays were detected from PSR B1957 -|- 20. 
It is out of the question that the pulsar is the only 
source of X-ray energy, but there are still lots of 
uncertainties in the explanation of the actual phys- 
ical mechanism. There arc no data so far on thc X- 
ray spectrum, and it is still even not clear whether 
these X-rays are emitted by the millisecond pulsar 
itself, or they originate due to physical processes 
taking place somewhere between the pulsar and 
the white dwarf. Thompson (1995) points to the 
fact that thc large spin-down luminosity of PSR 
B1957 -|- 20, taking into account the radius of the 
companion, gives the energy flux at the compan- 
ion's orbit which is seven times more than the one 
existing at the surface of our sun. In the model of 
Arons & Tavani (1993) the pulsar wind, shocked 
by its interaction with the wind off the compan- 
ion, emits X-rays and 7-rays with a soft X-ray 
luminosity comparable to that observed. A large 
fraction of this radiation is then absorbed by thc 
companion's surface facing the pulsar, due to the 
combination of relativistic beaming and the prox- 
imity of thc companion to thc shock. They show 
that the observed luminosity of X-rays can easily 
power also companion's optical emission. 

To the best of our knowledge, there is so far 
no spectroscopic data on the optical counterparts 
of the companion stars in both binaries. Let us 
first assume that the optical emission will be found 
to be nonthermal. The fastest "beam" particles 
of the pulsar wind with 7j_ ~ 10^~^ can reach 
the strong magnetic field at low altitudes from 
the stellar surface of the white dwarf, retaining 
their ultrarelativistic energies. Estimations show 
that such particles, gyrating around magnetic field 
lines, can emit a synchrotron emission falling into 
optical band. This might occur very close to the 
stellar surface, at distances from the companion's 
center that compare with the companion's opti- 
cal sizes calculated from the observations (see Sec- 
tion 1). 

Another mechanism can be suggested to explain 
the companion's high-energy radiation if the lat- 
ter proves to be of thermal nature. Relativistic 
electrons and positrons trapped by the magnetic 
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field of the white dwarf should "bounce" between 
higher-field points near the North and South poles 
(in the same manner as in the ordinary magnetic 
trap or in the Earth's magnetosphere) and slowly 
precess around the star due to V-B and curvature 
drift. The most energetic particles with the ratio 
v\\/v±_ at white dwarf's equatorial plane 



V± \ -Bmin 



1 



1/2 



(23) 



are in a so called "loss cone", and should be 
"poured out" on the magnetic poles. It can be 
demonstrated that broad polar regions of the com- 
panion star will be heated up to temperatures high 
enough to power the thermal X-ray emission. If 
this happens, the thermal optical emission can be 
just a long-wavelength tail of the Planck X-ray 
spectrum. 

We would like to notice that the same physi- 
cal mechanism was found to operate in the case 
of enigmatic pulsar Geminga (Gil et al. 1998). 
Namely, it appeared that the cyclotron damping 
of radio waves in the magnetosphere of this pulsar, 
combined with the almost aligned geometry, leads 
to the apparent absence of its emission at radio 
frequencies higher than about 100 MHz. More- 
over, we found that the cyclotron resonance has 
a strong impact on the general problem of the es- 
cape of radio waves from pulsar magnetospheres 
(Khechinashvili 1999). We plan to develop the 
eclipse model in forthcoming papers, where we ex- 
pect to achieve better agreement with the observa- 
tions and also include other eclipsing binary sys- 
tems into the framework of our model. 
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